Expansion velocity of a one-dimensional, two-component Fermi gas during the sudden 

expansion in the ballistic regime 
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We show that in the sudden expansion of a spin-balanced, two-component Fermi gas into an 
empty optical lattice induced by releasing particles from a trap, over a wide parameter regime, the 
radius R n of the particle cloud grows linearly in time. This allow us to define the expansion velocity 
V ex from R n — V e *t. The goal of this work is to clarify the dependence of the expansion velocity 
on the initial conditions which we establish from time-dependent density matrix renormalization 
group simulations, both for a box trap and a harmonic trap. As a prominent result, the presence of 
a Mott-insulating region leaves clear fingerprints in the expansion velocity. Our predictions can be 
verified in experiments with ultra-cold atoms. 
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I. INTRODUCTION 



Research into the non-equilibrium properties of 
strongly correlated many-body systems has emerged into 
a dynamic and active field, driven by the possibility to 
address questions such as thermalization [l|, 0] , the prop- 
erties of steady states, or state engineering in ultra-cold 
atomic gases While substantial theoretical attention 
has been devoted to quantum quenches in homogeneous 
systems 0, more recently, set-ups that give rise to fi- 
nite particle or spin currents have been studied as well, 
both from the theoretical side (I1-[TE| and in experiments 
(see, e.g., Refs. Il6l - [2ll ). Using these approaches allows 
one to investigate transport properties of strongly corre- 
lated many-body systems - in and out-of-equilibrium - in 
cold atomic gases that arc of great interest in condensed 
matter theory. 

Our work is motivated by the experiment by Schnei- 
der et al. [17f who have studied the expansion of a 



two-component Fermi gas in an optical lattice in two 
and three dimensions (described by the Fcrmi-Hubbard 
model [111, [23| ) , starting from an almost perfect band in- 
sulator. The qualitative interpretation of their results is 
that, besides a ballistically propagating halo of particles, 
at finite interaction strengths a core of diffusively expand- 
ing particles exists [l7|. In the case of one-dimensional 
(ID) bulk systems relevant for condensed matter prob- 
lems and on the level of linear response theory, ballistic 
dynamics of interacting particles can be traced back to 
the existence of non-trivial conservation laws 2A\ . For in- 
stance, the fact that the energy current is conserved for 
the ID Heisenberg model renders its spin transport ballis- 
tic away from zero total magnetization |24l42a |. whereas 
at zero magnetization there exists a quasi-local quantity 
[27j , which is conserved only for the infinite system, that 
gives rise to ballistic dynamics. While for the ID Hub- 



bard model, the understanding of its transport properties 
is by far less complete than for the Heisenberg chain, one 
might be tempted to expect similar quantities to play a 
role for the latter model as well (24| . 

A qualitative difference between the sudden expansion 
in an optical lattice compared to steady-state transport 
measurements in condensed matter systems is that, in 
the latter case, the background density determines trans- 
port coefficients, whereas in the former case, the density 
itself becomes time-dependent (l7| and all particles par- 
ticipate in the dynamics. As a consequence, in diffusive 
regimes, the dependence of the diffusion coefficient on 
density needs to be accounted for. In the ballistic case, 
as we shall see, the expansion velocity always depends on 
all momenta that are occupied in the initial state and not 
on just those close to the Fermi wave- vector. Therefore, 
a parameter regime complementary to condensed matter 
systems can be accessed with cold atoms. 

Theoretical results for the expansion of interacting 
bosons or fermions in optical lattices are mostly avail- 
able for the ID case, for which exact numerical methods 
give access to at least the short time dynamics via the 
adaptive time-dependent density matrix renormalization 
group (tDMRG) method [28l - l3ll ] or exact diagonaliza- 
tion (ED) H Hf. The richness of the non-equilibrium 
physics encountered in the expansion manifests itself in 
the observation of the dynamical emergence of coher- 
ence [1, H, [HI [Hj], which, for bosons, leads to the phe- 
nomenon of dynamical quasi-condensation [H, fill |32| and 
the intriguing phenomenon of the fermionization of the 
momentum distribution function (MDF) (El ITEl l33l |34| . 
In the case of a two-component Fermi gas, the short- 
time dynamics of the MDF and correlation functions [1] , 
the emergence of metastable states 0, and the time- 
evolution of density profiles for specific initial conditions 
have been investigated d, H, HEtHtJ ■ 

In the present work wc study the ID Hubbard model 
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FIG. 1: (Color online) Box trap: Typical contour plot of the 
density (rii(t)) during the expansion from a MI (U = 8 J, 
"-init = 1, iinit — 20). The slanted lines indicate the speed 2 J 
at which the MI melts. 




FIG. 2: (Color online) Box trap: Radius R n (t) for initial 
densities n ini t = 0.2, 0.4, 0.8, 1 at U = 8J and L init = 40 
(corresponding to N = 8, 16, 32, 40). 



and we concentrate on the sudden expansion starting 
from initial states that are Mott insulators (MI), i.e., that 
have an integer filling of n- m n = 1, Tomonaga-Luttingcr 
(TL) liquids (ninit < 1), or systems in a harmonic trap. 
In the latter case, depending on filling and interaction 
strength, several phases may coexist in separate shells 
[38j . We analyze the dependence of the expanding cloud's 
radius R n (t) on time t and search for conditions to obtain 
ballistic dynamics, for which R n {i) = V cx t is a necessary 
criterion. In that case, the expansion velocity I4 X is a 
well-defined quantity, and, as a key result of our work, 
we clarify its dependence on the initial conditions. 

Our main results are: (i) In the regime of low densities, 
i.e., njnit < 1) we observe a linear growth of the cloud's 
radius with time, allowing us to define V ex . (ii) In general, 
the expansion speed V ex depends in a non-monotonic way 
on the initial density. In the case of the expansion from 
a MI, V^j X = y/2J, independently of U. (iii) Our findings 
are robust against the presence of a harmonic trap in the 
initial state. 

Note that, in a generic system, one expects ballistic 
dynamics in the long time limit, where the gas becomes 
so dilute that interactions cease to matter. Here we show 
that ballistic dynamics sets in immediately after the gas 
is released from the trap when the density is actually still 
comparable to the initial density. 

The structure of the paper is the following: In Sec- 
tion nn we introduce the model and define the radius of 
the cloud. Section ITO! discusses the expansion from a box 
trap, i.e., starting from a homogeneous density. We first 
show that the dynamics is ballistic by analyzing the ra- 
dius and the particle currents and second, we present a 
detailed analysis of the expansion velocity as a function 
of density and interaction strength. In Sec. IIVI we test 
our findings against the inhomogeneity introduced by a 
harmonic trap. We summarize our findings in Sec. [V] 
In Appendix A, we discuss the diffusion equation in one 
dimension. Appendix B contains a finite-size analysis of 
the expansion velocity for various cases. 



II. MODEL AND SETUP 



Our study is carried out for the ID Hubbard model: 
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c\ a is a fermionic creation operator with spin a =t, 4- ac t~ 

ing on site i, rii a — c\ a Ci a , rii = ^Za nia - U 1S ^ ne ons ite 
repulsion, and J, is the hopping matrix element. Open 
boundary conditions are imposed, L > 100. is the num- 
ber of lattice sites, and N the number of particles. We set 
h and the lattice spacing to unity and thus measure time, 
velocity and particle current in the appropriate units in 
terms of the hopping matrix element. 

We prepare initial states as the ground state of H = 
Hq + H con f Q. We consider two cases: First, the expan- 
sion from a box trap (i.e., (rij) ^ for i± < % < i 2 ; 
(i 2 - i-i) = Ljnit, ninit = A/Linit) enforced by using 
H con { = J2i £ i n i with a large ej > U for i < ii;ia < i and 
zero otherwise). The second example is the expansion 
from a harmonic trap, for which H con { = V^2 i {i — io) 2 rn. 
We turn off iJ con f at t = 0. In our tDMRG runs, we use 
a Krylov-space based method [3^, |4(| , with time steps of 
St J = 0.25 and we enforce a discarded weight of 10~ 4 or 
smaller. 

The main quantity of interest is the radius of the par- 
ticle cloud that we define via 
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if(n)(>-j,P-l5(i = G}. (2) 



For the expansion from a box, io = L/2 + 0.5. 
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III. EXPANSION FROM A BOX TRAP 

We first discuss this idealized case to avoid the com- 
plication of dealing with particles originating from dif- 
ferent shells, as would be the case with a harmonic trap 
(note, though, that box-like traps can also be generated 
in experiments jU |42j ) ■ A typical example for the time- 
evolution of the density (m(t)) is shown in Fig. [T] for 
the expansion from a MI with U = 8 J. The MI melts 
on a time scale of i me it 'S Anit/(2J), where 2J is the 
largest possible velocity in the empty lattice, since the 
single-particle dispersion is tk = — 2Jcos(fc) For 
t > t m eit-, two particle clouds form that propagate into 
opposite directions, yisible as two intense jets (compare 
Refs. 0,0, III Eli!]). 

In Fig. [2 we display the radius R n (t) at U = 8 J for 
various initial densities at U = SJ. Clearly, for < 1, 
Rn(t) = V CK t. We stress that R n {t) ~ t sets in imme- 
diately after the gas is released from the trap. This in- 
cludes, in particular, the expansion from a MI at any U, 
while for ninit > 1, the radius deviates from R n (t) ~ t 
0. Based on the observation of R n (t) ~ t on short and 
intermediate times, when local densities arc still large, 
together with the fact that interacting particles behave 
similar to non-interacting ones (which, in the absence of 
disorder, expand with R n ~ t), we classify the dynamics 
as ballistic. 

In our situation, the notion of ballistic dynamics is 
strongly corroborated by analyzing the time dependence 
of the total particle current in each half of the system, 

Jl/2 = Ei>z,/2^ = - iJ Y.A c \+i<j c ™ ~ h.c.)], which 
is shown for U/J= 2,8 in Fig. [3] After the two jets in 
Fig. Q] are well separated from each other, Jl/2 takes a 
constant value, which we consider a hallmark feature of 
ballistic dynamics [Hj]. 

However, in one dimension, there is a subtlety as cer- 
tain solutions of the diffusion equation can also give rise 
to a linear increase of the radius with time (if properly 
defined). Such a scenario happens in the dilute limit 
(which we do not study here), yet it results in a strong 
dependence of the expansion velocity on the total parti- 
cle, which is clearly different from our case as we shall 
see below. Further details are given in Appendix A. 

The observation of a linear increase of the cloud radius 
with time implies that V ex should be fully determined by 
properties of the initial state, such as the MDF, energy 
per particle, or density. In the non-interacting case, this 
is obvious, since V ex can be calculated from the knowl- 
edge of the MDF. To guide the interpretation of the in- 
teracting case and to understand the dependence of V ex 
on U and n; n ; t , we next study the two exactly solvable 
limits U = and U = oo. 



A. Box trap, V c * at U — 

At U = 0, opening the trap simply means that particles 
will propagate with a velocity Vk = 2 J sin(fc) with a prob- 
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FIG. 3: (Color online) Box trap: Total particle current in 
each half of the system as a function of time, i.e., Ji/2(t) := 

Si>i/2 Ji ' for ( a ) U = 8J ' ( b ) U = 2J and n ini t = 
0.2,0.4,0.6,0.8,1. We observe that after some transient dy- 
namics, (Jl/2(£)} = const, supporting the observation of bal- 
listic transport. 

ability given by the MDF in the initial state, which 
is n k = {l/N)Y jl ^ <T e-^- m ^ k {c\ <jCma ) . The momenta 
are chosen to match the open boundary conditions in the 
box, i.e., k = L ; / = l,...,L init . By a straight- 
forward evaluation of R„(t) from Eq. ([2]) and using the 
time-dependence of creation and annihilation operators, 
known exactly at U = 0, we obtain as the average 
velocity of all particles in the initial state: 

K 2 x = ^£^n fe . (3) 

A; 

In the [7 = case, the initial MDF thus completely de- 
termines the expansion velocity. However, this is an over- 
complete set of constraints: For a very large N, where 
boundary conditions cease to matter, wc can evaluate 
Eq. ([3]) analytically: 

K 2 X = 2J 2 [k F - cos(fc F ) sm(k F )]/k F , (4) 

which yields the full dependence on the initial density 
at U = through kp oc Hinit alone. We can interpret 
Eq. (U]) in two ways: If U = 0, k F = 7rrii n it/2, whereas 
for U = oo, kp = 7rni n i t . Using ED, we have verified 
the validity of Eq. (UJ by extracting V ox from the time- 
dependence of R n (t) for N ~ 160 (see Fig.|8]in Appendix 
B). 

B. Box trap, V cx at U 

In the interacting case, we extract the expansion ve- 
locity V^x from the tDMRG data (i.e., the slope of curves 
such as the ones shown in Fig. [2]). The results for selected 
values of U arc collected in the main panel of Fig. @] 
(symbols). Wc emphasize four main observations: (i) 
For the expansion from the MI, we obtain Ve X = V%J at 
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FIG. 4: (Color online) Box trap: Main panel: V" cx vs ?ij n it at 
U/J = 0,0.5,2,8,oo for Linit = 20 [see the legend, symbols 
are tDMRG, the solid and the dashed lines are derived from 
Eq. ©]. Inset: V cx from tDMRG (open squares: U = 8J, 
open triangles: U = 2J) vs. Kef (solid symbols) from 
non-interacting reference systems at a finite temperature (see 
Sec.|niC]for details) 



any U > 0. (ii) At a fixed density, V cyi increases mono- 
tonically with U. (iii) For U > 4J, the maximum of 
the expansion velocity is at an incommensurate density 
0.5 < riinit < 1. (iv) The expansion velocity is always 
very different from characteristic velocities of the initial 
state and much smaller than 2 J, the largest possible ve- 
locity. It is also much smaller than the charge velocity 
pfjj at small densities and at n; n it = 1, where the charge 
velocity drops to zero, V cx remains finite. 

At U = 0, the first observation is a consequence 
of particle-hole symmetry, reflected in the MDF: n k 
is point-symmetric about the point (kp = -k /2,nk F ). 
Since vl F+Skp = vl F _ 5kp , from Eq. ©, we conclude 
V cx = \/2J. The MDF at U > has the same symmetry 
property, hence we expect a similar behavior, confirmed 
by tDMRG. Of course, Eq. does not directly apply to 
the interacting case. Since the total energy Ejj = (Ho) 
is conserved, for U > 0, Eq. © is incompatible with this 
initial condition set by U > and roi n it. However, we 
shall see that the observation of V ex = \2J for U > 
can also be understood as a consequence of symmetry 
properties. 

We can further use the exact result Eq. (|H) to explain 
the observations (ii)-(iv). The U = and the U = oo re- 
sult are the solid and the dashed lines in the main panel 
of Fig. |4l respectively, and therefore, increasing U from 
U = to U = oo at a fixed density simply takes us from 
the limit of a non-interacting two-component Fermi gas to 
the limit of non-interacting spinless fermions. To under- 
stand that the maximum of V eyi is at an incommensurate 
?iinit for U > 4 J, one needs to take into account that on 
the one hand, in a ID cosine band, the maximum veloc- 
ity is at k = 7r/2, but on the other hand, the density of 
states takes its minimum there. As a consequence of this 
competition, i.e., the decrease of vs the increase of the 
density of states as one moves away from ir/2, the largest 
expansion velocity is at ninit !■ Finally, property (iv) 
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FIG. 5: (Color online) These plots illustrate how the non- 
interacting reference systems are constructed using Eqs. (5), 
(6) and (7) for the example of ninit = 0.4 for (a),(b),(c): 
iinit = 40 and (d): Linit = 160: (a) Temperature dependence 
of the chemical potential at fixed ninit = 0.4. (b) Total energy 
from Eq. (6) as a function of temperature, (c) V^f as a func- 
tion of temperature. The horizontal lines in (b) are DMRG 
results for E\j for the initial states used in the expansion at 
zero temperature for (7 = 0, 2 J, 8 J, while in (c) they mark the 
resulting V T 2 ci = V T 2 ct (E,n). (d) As T increases, V r 2 ci 2J 2 
from below, i.e., this is the largest expansion velocity that the 
reference systems can produce for n; n it < 1. 



is a consequence of all particles propagating and not just 
those with momenta close to kp. 

On a technical note, we have checked the dependence 
of Vox on particle number, keeping ni n ; t = N/L- m i t fixed. 
Finite-size effects are the largest at small initial densities, 
yet for densities rii n i t > 0.5, our tDMRG results obtained 
with Ljnit = 40 show little quantitative differences com- 
pared to smaller L;„it and V cx becomes independent of N 
as shown in Appendix B. 



C. Reference systems 



It is now a compelling question to ask how many con- 
straints suffice to determine the expansion velocity From 
the solution of the non-interacting case, we conclude that 
density and energy are relevant quantities. To check this 
conjecture for the interacting case, we construct non- 
interacting reference systems that are at a finite tem- 
perature 47J. The temperature is chosen such that the 
reference system has the same energy as the interacting 
system and the same particle number, and it lives in the 
same box potential of length -Lhiit. 
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Hence we solve this set of equations: 

N = ^/(e fe -/i,T), (5) 

A:. (J 

Eu = £> fe /(e fc -M,T), (6) 

A:. (J 

where f{x,T) is the Fermi function. We proceed as illus- 
trated in Fig. [5] For a given U and AT, we compute the 
total energy Eu in the initial state with DMRG. First, 
we find the chemical potential p, = p(T) from Eq. ([5|. 
which only depends on N. Using this p(T) curve, we 
determine the pair of (/z,T), for which we get the right 
energy Ejj. From these results, Eq. yields the expan- 
sion velocity V re [ of the reference system. Obviously, the 
maximum velocity that these reference systems, which 
have the dispersion = — 2Jcos(fc) of the empty lattice, 
can produce is V lc f = y/2J at any density n; n ;t < 1 as 
T — > oo. Within that constraint, the agreement between 
V ex and our reference systems is excellent, as we illustrate 
for U/J = 2 and 8 in the inset of Fig.|4j Apart from those 
densities for which, at U = 8J, V cx > V2J, V le { m V cx 
within our numerical accuracy. In the particular case of 
"init = 1, our reference systems also yield V Te { = s/2J in- 
dependently of Eu, consistent with the tDMRG results 
of Fig. SI This is a consequence of the aforementioned 
symmetry property of the MDF, which also applies to 
T > 0. 



IV. EXPANSION FROM A HARMONIC TRAP 

Our results so far establish a relation between proper- 
ties of the initial state and the expansion velocity that 
could be probed in experiments. We next test the robust- 
ness of our predictions for V ex = V ex (U, n; n it) against the 
inhomogeneity induced by a harmonic potential. 

We focus on three types of initial states: (i) Only a 
TL, i.e., (rii) < 1 in the entire trap, (ii) a MI shell in the 
center, surrounded by TL wings, and (iii) a three-shell 
structure with an incommensurate density in the center 
> 1, surrounded by first, a MI shell and second, a 
TL shell with (rij) < 1. For a given U > 0, these regimes 
are separated by critical characteristic densities p\ and 
p2, where p = Ny/V/J is the effective density in a system 
with a harmonic trap [H, HH . 

For all three cases, R n (t) is shown in Fig.[B]for U/J = 2 
and 8. We observe that, after releasing the particles from 
the harmonic trap, the cloud still expands with R n {t) ~ t 
in cases (i) and (ii), i.e., R n (t) ~ t [see Fig.|6](a) and (b)] 
whereas in case (iii), the increase of the radius is slower 
than linear in t [see Fig. [6] (c) and (d)]. In that regime 
and for U > 4 J, the system can be viewed as a mixture of 
single atoms propagating with velocities Vk ~ J and two 
fermions repulsively bound into a doublon, which, due to 
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FIG. 6: (Color online) Harmonic trap: Radius of the particle 
cloud for the expansion from a harmonic trap for (a) N = 
10, 30 at U = 8 J, (b) N = 10, 20 at U = 2J, (c) TV = 48 
at U = 8J and (d) N = 30 at U = 23. For U = 8J, these 
parameters correspond to the initial states shown in the insets 
of Fig. [7] The solid lines are fits to the data, dotted lines are 
guides to the eye. We find R n (t) ^ t whenever densities in the 
center of the trap are larger than one and R n (t) ~ t otherwise. 

energy conservation, does not decay on time scales cx 1/J 
and is much slower with typical velocities ~ J 2 /U 
[HI . For illustration, the values of p\ and p2 as well as 
typical density profiles are included in Fig. [7] for U = 8J 
(vertical lines and lower insets, respectively). As is evi- 
dent from Fig. [71 the overall dependence of V cx = V ex (p) 
resembles that of the expansion from a box trap, with 
a maximum in V cx emerging as U > 4J. Most impor- 
tantly, as soon as the MI forms in the center of the trap, 
indicated by the vertical solid line at p = p\, the expan- 
sion velocity approaches a constant value at V cx > ^/2J 
from above. The contribution to V^x of low-density shells 
surrounding the MI is suppressed by increasing U or p 
since both favor a large relative fraction of all particles 
in the MI shell to minimize the contribution from the in- 
teraction energy. In contrast to the expansion from the 
box, the limit of U — oo (dashed line) is approached very 
slowly since the shell structure in a trap depends strongly 
on U and p. 



V. SUMMARY 

We studied the sudden expansion of a spin-balanced 
two-component gas in ID, released from a trap. Our 
main results are two-fold: First, the cloud expands bal- 
listically as long as initial densities are small, including, 
in particular, the MI state. Second, the expansion ve- 
locity, defined through R n (t) = V cx t strongly depends on 
initial density and thus, its measurement can provide in- 
formation on the initial state. For instance, deviations 
from our predictions could indicate the presence of de- 
fects in the initial state preparations. Our quantitative 
predictions can be tested in an experiment that realizes 
the set-up of Ref. [l7| in ID. 

Furthermore it would be interesting to study the ra- 
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FIG. 7: (Color online) Harmonic trap: V ox vs p (U/J = 
2,8,40,oo, V = 0.016J). The vertical solid line marks the 
formation of a MI shell in the trapped system at pi and the 
vertical dashed line the formation of a core with (m) > 1 at 
p2, both for U — 8 J. The horizontal line is V cx = y/2J. Sym- 
bols are tDMRG results, the dashed line was obtained from 
ED. We have verified that the results are remarkably stable 
against changes in the particle number at fixed p by produc- 
ing the V C x = Kx(p) curve at a different V (see Fig. 10 in 
Appendix B). Lower insets: typical initial density profiles in 
the regimes p < pi, pi < p < p2 and p2 < p for U = 8 J. 



dius of an expanding cloud and the expansion velocity for 
other experimentally relevant systems such as the Bosc- 
Hubbard model or spin imbalanccd mixtures. While we 
have presented phenomcnological evidence for ballistic 
dynamics, we have here not touched upon a potential 
relation with integrability and non-trivial conservation 
laws [24[ , leaving this for future research. It also remains 
as an open question to identify interacting models in one 
dimension and parameter regimes in which diffusive dy- 
namics dominates during the sudden expansion, which 
might be challenging since even non-integrable models 
may have very large conductivities (see, e.g., Ref. [26j). 
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Appendix A: Linear increase of the radius from a 
nonlinear diffusion equation 

Here we discuss solutions of the diffusion equation in 
one dimension in the limit of a very dilute gas. Since the 
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init 



FIG. 8: (Color online) Box trap: Comparison between the 
exact result for V ax [Eq. (4) of the main text, valid in the 
limit of large particle numbers] and exact diagonalization in 
the limits U — and U — oo. ED data are shown for £mit = 
40 and Linit = 160. In the latter case, we find very good 
agreement between the analytical result (lines) and the ED 
data (symbols) for all densities. 



sudden expansion scenario considered in this paper in- 
volves the propagation of all particles, the dependence of 
the diffusion constant D on the local density n(x,t) be- 
comes relevant, and as a consequence, the relevant diffu- 
sion equation is in general a nonlinear one (see, e.g., [l7|). 
Focussing on the very dilute limit we use D ~ l/n(x,t) 
(see the discussion in Ref. [13, El))- The resulting diffu- 
sion equation (after rescaling of the time variable ): 



d t n(x,t) = d x -d x n. 
n 



(8) 



has a self-similar solution with particle number conser- 
vation in ID El: 



n(x, t) 



2t 



c 2 + v 2 t 2 



(9) 



First of all, one realizes that our definition of the ra- 
dius R n (t), Eq. ©, cannot be used here. In the analysis 
of experimental data, it is common practice to define the 
radius as the half-width at half-maximum of the expand- 
ing cloud (l7|. Using this definition, the solution Eq. © 
yields indeed R n {i) = vt, similar to the ballistic dynamics 
discussed in our work. We would like to stress, though, 
that the sudden expansion described in the main text 
is genuinely different in some important respects. First, 
Eq. ijSJ) is only valid in the dilute limit while the time- 
dependent DMRG gives us access to short and interme- 
diate time-scales only where the gas is not necessarily a 
dilute one yet. Second, Eq. ((9]) is a solution for which 
the expansion velocity v depends strongly on the parti- 
cle number via TV = J_ n(x, t)dx = 2n/v, which is not 
observed in our case (compare Fig. [4] and |9]). Based on 
these differences, we conclude that diffusive dynamics is 
very unlikely to be realized for the ID Hubbard model in 
the sudden expansion. 
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FIG. 9: (Color online) Box trap: Expansion velocity as a 
function of initial density for U — 8J and different iinit 
I i ) . 20, 40. For rii n it > 0.6, finite-size effects are remarkably 
small. 
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FIG. 10: (Color online) Harmonic trap: Expansion velocity 
for V = 0.008J and V — 0.016J as a function of effective 
density p = N ■ y/V/J for U = 8J. The expansion velocity 
is remarkably stable against changing the particle number at 
fixed p. 



Here wc address the question of how our results for the 
expansion velocity depend on the overall particle number 
at a fixed density nj n it- First, we consider the box trap 
and we compare our analytical result for large TV [Eq. 
to exact diagonalization in the nonintcracting limits in 
Fig. [SJ For TV = 40, we find good qualitative agree- 
ment with small finite-size effects, which are the most 
pronounced for riinit < 0.5. For TV = 160, the deviations 
between the analytical expression for TV — > oo and data 
for a finite TV are already barely visible except for very 
low densities. Second, we study the interacting system 
expanding from different box traps with Lj n i t = 10, 20, 40 
at a fixed density for U = 8 J. Figure [5] shows V cyi as a 
function of density. As in the non-interacting case finite- 
size effects are remarkably small whenever nj n i t > 0.6 
even for the smaller particle numbers. Finally, we turn 
to the expansion from a harmonic trap and analyze V ex 
for two different trapping potentials, V = 0.008J and 
V = 0.016 J. Fi g. [TOl shows V eyi as a function of effective 
density p = TV \/V/ J. Wc find that the expansion veloc- 
ity is very robust against changing the particle number at 
a fixed p. Overall, our results for the expansion velocity 
exhibit only minor finite-size effects in all studied cases. 
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